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Abstract 

We present efficient algorithms for detecting central and mirror symmetry for the 
case of algebraic curves defined by means of polynomial parametrizations. The al- 
gorithms are based on an algebraic relationship between proper parametrizations 
of a same curve, which leads to a triangular polynomial system that can be solved 
in a very fast way; in particular, curves parametrized by polynomials of serious 
degrees/coefficients can be analyzed in a few seconds. In our analysis we provide a 
good number of theoretical results on symmetries of polynomial curves, algorithms 
for detecting rotation and mirror symmetry, and closed formulae to determine the 
symmetry center and the symmetry axis, when they exist. Some observations and 
empiric results for the case of polynomial parametrizations with floating point co- 
efficients are also reported. 



1 Introduction 



This paper deals with the problem of detecting the symmetries, if any, of a 
curve denned by means of a polynomial parametrization, i.e. a pair (x(t), y(t)) 
where both x(t) and y(t) are polynomials. The question of finding the sym- 
metries of an algebraic curve has been previously investigated mainly because 
of its applications in pose estimation and patter recognition. Essentially, the 
problem is the following: a situation that is often studied in pattern recognition 
is how to choose, from a database of algebraic curves representing different 
objects, the one that best fits a given object, also represented by an alge- 
braic equation; this question is addressed, among many others, in [I] (where 
the database simulates different aircraft prototypes), [7], [10], [13] (where the 
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database corresponds to silhouettes of sea animals), [H] or [15]. In order to do 
this one has to detect if the curve to be recognized is in fact the result of ap- 
plying an affine transformation (typically translations, rotations, etc.) to some 
curve in the database; in turn, this implies to set the curve in a "canonical 
position" that makes recognition possible. The problem of detecting symme- 
tries comes along then as a means for setting properly the curve. In particular, 
this question has been addressed, for example, in [4j, using splines, in [I], [2], 
[T7j . by means of differential invariants, in [5], [BJ, [32], using a complex repre- 
sentation of the implicit equation of the curve, or in jij, [TT], p2], [16], using 
moments. 

It is also worth observing that in most of these papers it is assumed that the 
curve is given by means of its implicit equation. Exceptions to this are [3J, 
where the curve is assumed to be a spline (i.e. a union of pieces of polyno- 
mial parametrizations), and pQ, [2], and other several papers on differential 
invariants, where the input is considered to be a parametrization without any 
restriction on its functional form. Furthermore, in almost all the papers on 
the question, the algorithm that is provided to detect symmetries is basically 
numerical, and therefore the output is approximate, and not exact. This is 
not a problem when the form to be recognized is, up to a certain extent, 
fuzzy, or when there are missing data, common situations in pattern recog- 
nition. But it is undesirable if the input is exact. Up to our knowledge, the 
exceptions to this are the papers on differential invariants ([TJ, [2], etc.), and 
[5], [BJ. However, when applied to produce an exact output, the former can 
only deal with parametrizations of low degree (in fact, the ultimate idea in 
the papers on differential invariants is to adapt the underlying theory to a 
numerical framework). On the other hand, the latter, i.e. [5], [BJ, do provide 
deterministic algorithms for implicit algebraic curves f(x,y) = 0, and yield 
exact answers in an efficient and elegant manner (even for serious curves). 
The only limitation is that these algorithms require the absolute values of the 
leading coefficients with respect to the variables x, y of f(x, y) to be different, 
which is a limitation in certain cases. 

In this paper, we address the problem from a different perspective. On one 
hand, we focus on curves defined by means of polynomial parametrizations. 
Even if this seems too restrictive, polynomial parametrizations are widely 
used, for example, in Solid Modeling and Computer Aided Geometric Design, 
and stay at the core of the notion of spline curve, commonly used in many 
applications. Furthermore, the algorithms that we provide can be directly 
conducted from the parametrization, and therefore do not require to compute 
the implicit equation of the curve (costly or even impossible for high degrees); 
also, no condition on the leading coefficients of the curve is needed. Finally, 
we assume that our input is exact (the coefficients of the parametrization 
are required to be real numbers, not necessarily rational) and we provide 
also exact, i.e. deterministic, algorithms for checking whether the curve is 
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symmetric or not, and for determining the elements of the symmetry, in the 
affirmative case. 

The algorithms that we provide are very efficient and are able of analyzing 
curves with serious coefficients and/or high degrees in just a few seconds. 
The idea behind these algorithms comes from Computational Real Algebraic 
Geometry, and exploits an algebraic property linking two "good" (in a sense 
that is introduced in Section [2j parametrizations of a curve. By applying 
this property, the analysis of symmetries is reduced to solving polynomial 
systems which are triangular, and that can therefore be analyzed in a fast and 
efficient way. Additionally, the development of our results produce a number 
of theoretical results on the symmetries of polynomially parametrized curves, 
and closed formulae for the symmetry center and the symmetry axis, in the 
cases when they exist. 

So, in the sequel we will analyze rotation and central symmetry (i.e. symmetry 
with respect to a point) in Section[2j and mirror symmetry (i.e. symmetry with 
respect to a line) in Section [3] In both cases we report examples and timings 
showing the efficiency of our algorithms. In Section [4] we report some observa- 
tions and empiric results on polynomials with floating point coefficients, i.e. 
coefficients known up to a limited precision. We finish the paper with a section 
on conclusions and further work. 



2 Rotation and Central Symmetry 

Along the paper, we let C be an algebraic curve admitting a real polynomial 
parametrization (p(t) = (x(t),y(t)), i.e. x(t),y(t) e R[t] (we will summarize 
this by saying that C is a real polynomial curve). Furthermore, we will also 
require ip(t) to be proper, i.e. injective for almost all values of t. For example, 
{t,t 2 ) is a proper parametrization of a parabola, while (t 2 ,t 4 ) is not (notice 
that whenever t moves over the complex numbers, the latter parametrizes the 
parabola as well, and almost all the points of the parabola are generated by two 
different values of t). The properness of a parametrization is considered to be a 
good property, since non-proper parametrizations provide a "redundant" rep- 
resentation of the curve, that might cause problems when plotting the curve, 
intersecting it with other curves, or simply determining notable points of it. 
From Theorem 6.11 in [0], it is guaranteed that a polynomial curve can always 
be properly and polynomially reparametrized; more than that, we can ensure 
that in fact one can always find a proper polynomial parametrization over the 
reals. Indeed, if (p(t) = (x(t),y(t)) is polynomial but it is not proper, from the 
algorithm in page 193 of [9] we have that it can be properly reparametrized 
without extending the ground field; furthermore, if this reparametrization is 
not polynomial, then the algorithm in page 199 of [9] leads to a polynomial 
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reparametrization, again without extending the ground field. So, in the rest 
of the paper, without loss of generality we will assume that ip(t) is proper. 

2.1 Rotation Symmetry and Prohibitions. 

We say that C has rotation symmetry iff there exists a point Pq G R 2 and an 
angle G (0, 2ir) such that the rotation of center P and angle leaves C 
invariant. In the special case when <fi = ir, we say that C has central sym- 
metry, i.e. that C is symmetric with respect to Pq (the center of symmetry). 
The following theorem proves that in our case, the only form of rotation sym- 
metry that a polynomial curve can exhibit is central symmetry. Here we use 
the notion of infinite branch of a parametric curve: by this notion, we mean 
intuitively a part of C that goes to infinity; for example, the parabola (t, t 2 ) 
has, according to this terminology, two infinite branches, one for t = oo and 
another one for t = — oo. In fact, every polynomial curve has only two infinite 
branches; however, rational curves may have more than two (because there 
may be t- values where the denominator of either x(t) or y(t) vanishes). 

Theorem 1 Real polynomial curves cannot have any other form of rotation 
symmetry other than central symmetry. 

Proof. If C would present rotation symmetry other than central symmetry 
then it would consist of the union of m different copies, with m > 3. Since 
each copy would provide at least one infinite branch, we would get at least m 
infinite branches. But since m > 3, and the number of infinite branches of C 
is 2, we reach a contradiction. n 

Furthermore, the following result provides a necessary condition for C to have 
central symmetry. 

Theorem 2 Let C be a real polynomial curve, and let ip(t) = (x(t),y(t)) be a 
proper, real, polynomial parametrization of C If C has central symmetry, then 
limt^^xit) and Umt^-ooX^) , /zr?2 t _^ 00 |/(t) and limt^-ooy^) , exhibit different 
signs. In particular, deg t (x(t)) and deg t (y(t)) must be both odd. 

Proof. If C has central symmetry then it must exhibit two different infinite 
branches, that go to infinity in opposite quadrants. The rest follows easily. □ 

2.2 Detecting Central Symmetry 

According to the preceding subsection, central symmetry is the only form of 
rotation symmetry that C can exhibit. So, let us address here the problem of 
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detecting this kind of symmetry in an efficient way. For this purpose, previ- 
ously we need the following result, that is crucial for us. 

Lemma 3 Let <fi(t) and <f2{t) be two polynomial and proper parametrizations 
of the same curve. Then there exists a linear function L(t) = at + (3, with 
a, (3 G K, such that fi(t) = y9 2 (at + P) ^nd q^O. 

Proof. By Lemma 4.17 in [9], (fi(t) and (f2{t) are both proper parametriza- 
tions of C if and only if each one can be obtained from the other by ap- 
plying a rational reparametrization = — ^ , where a, (3,^,6 G C and 

jt + o 

a5 — (3 , -f 7^ 0. However, since both ipi(t) and (pzif) are polynomial, from Lemma 
6.8 in [9] we deduce that the reparametrization £(t) is in fact linear. Let us see 
that a, (3 61. Indeed, if a, {3 G C then since (pi(t) = y3 2 (at + (3) and fi(t) is a 
real polynomial parametrization, we have that v 9 2(« + t(3) = y3 2 (at + /?)■ But 
then v?2(£) cannot be proper, because for every (real or complex) value of t, it 
holds that at + (3 and at + (3 generate the same point. Since <f2{t) is proper 
by hypothesis, then a, (3 G K. Finally, since in this case 5 = 1 and 7 = 0, then 
a<5 — /?7 = a ^ 0. □ 

Along the paper, it will be more convenient to write (f(t) in complex form as 



z {t) = x{t) + i • y(t) 



Now C has central symmetry, with center of symmetry z$ (zq corresponds 
to a point Pq = (xo,yo), written in complex form), if and only if for every 
t-value, z*(t) = —z(t) + 2z is also a point of C (see Figure 1), i.e. if and 
only if z*(t) is another parametrization of C (in complex form). Now one 
may check that z*(t) also corresponds to a real polynomial parametrization; 
furthermore, since z*(t) is the result of composing cp{t) with a bijective planar 
transformation (the symmetry with respect to Zo)> then it is also a proper 
parametrization. So, C has central symmetry if and only if z(t) and z*(t) are 
two proper parametrizations of the same curve. These observations, together 
with Lemma [3j give rise to the following theorem. 

Theorem 4 The curve C has central symmetry if and only if there exist Zq G C 
and a linear transformation L(t) = at + (3, with a, (3 G R, a ^ 0, such that 
z*(t)=z(L(t)). 

Now writing z{t) = c„t n + c n _it n_1 + • — h c , where q G C for % = 0, 1, . . . , n, 
and L{t) = at + (3, with a, (3 G M, the statement in Theorem El leads to the 
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z*(t) = z + [-(z(t) - z 
= 2zq — z(t 




Fig. 1. Central Symmetry 
following (n + 1) x (n + 1) polynomial system, that we denote as <S: 



-(co - zo) + z = c + ct/3 + c 2 (3 2 H h 



n-1 



+ 



-Cl 

-c 2 



Ci« + c 2 • 2a/3 H h c n _i(n — l)a/3 n 2 + c n na/3 n 1 

c 2 a 2 +■■■+ Cr ^-^^ a 2 ^ + c n ^±V/r- 2 



"C n -1 



C n _i« 



n-1 



+ c n na n x /3 



Since by Theorem [2] n must be odd and the case n — 1 is trivial (C is a 
line), in the rest of the section we will assume that n > 3. Now we will 
refer to the above equations as (0), (1), . . . , (n), respectively (i.e. the equation 
(k) is the one containing a k ). Observe that Theorem [4] is equivalent to the 
existence of a solution (z , «, 0), with of the system 5. This system 

has a triangular structure that can be exploited for solving the problem in an 
efficient way: first of all, notice that from the last equation we have a n = — 1, 
with n odd. So, the only possible real value of a is a = — 1. Additionally, from 
the equation (n — 1) we deduce that /3 = ~ 2c " -1 . Hence, if C is symmetric then 
this expression must yield a real number. So, the only possible value for (a, 0) 
is (—1, ~ 2 ^~ 1 ). Finally, since the coefficient of z in the equation (0) is —2, 
whenever (—1, ~ 2c "~ 1 ) fulfill the equations (1), . . . , (n) we can solve for zq in 
the equation (0), and obtain the symmetry center. We summarize these ideas 
in the following theorem, where a closed expression for the symmetry center 
(in the case when it exists) is provided. 

Theorem 5 The curve C has central symmetry if and only if (3 = ~ 2c ^~ 1 G R, 
and a = —l,/3 = ~ 2 i c "~ 1 satisfy the equations (l),...,(n) of the system S. 
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Moreover, in that case the center of symmetry is the point zq (in complex 
form) given by 

zo = c + ~( Cl p + c 2 f3 2 + ■■■ + c n f3 n ) 

Remark 1 From Theorem 5. 3 in JEjj, it follows that C can have at most just 
one center of symmetry. In our case, the existence of at most one symmetry 
center follows also from the above theorem. 

This theorem provides in a natural way an algorithm for checking central 
symmetry, and for determining the center of symmetry in the affirmative case. 
This is illustrated by the following example. 

Example 1. Consider the curve C parametrized by (x(t),y(t)), where 

x(t) = 2 + 2(2t + l) 23 - (2t + l) 13 + 2(2t + l) 11 + 2(2t + l) 5 - (2t + l) 3 + 2t 

< 

y(t) = -2(2t + l) 23 + (2t + l) 13 - 2(2* + l) 11 + 2(2* + l) 5 - (2t + l) 3 + 2t 

One can check that the curve is proper by applying Theorem 4-30 in J2J/. In 
fact, it is perhaps instructive to mention that we have constructed this curve 
by considering first ip(t) = (2t 23 — t 13 + 2t u , 2t 5 — t 3 + t), which is obviously 
symmetric with respect to the origin, then substituting t := 2t + 1, and finally 
applying the change of coordinates {X = 1+x+y, Y = — 1— x+y}, which is the 
composition of a rotation and a translation of vector (—1, 1); as a consequence, 
we obtain a curve which is symmetric with respect to the point (—1,1). Now 
one can check that the implicit equation f(x,y) of C has degree 23, that the 
infinity norm of this implicit equation is close to 2 600 ; and that the leading 
coefficients of f w.r.t. x and y are equal. This last circumstance makes it 
not possible to apply the method in |5]/; nevertheless, the size of the implicit 
equation suggests that it would be difficult to work in implicit form, anyway. 
Using our method and Maple 14, running in a computer with 8 Gb of RAM 
and a CPU revving up to 2 GHz., it takes 0.374 seconds to construct the system 
S and check that a = —1, (3 = —1, zq = 1 — i is a solution. So, we recover 
(1, —1) as the symmetry center of the curve. 

In the case when c n _i = the analysis of the system S is easier. Indeed, in 
this case from the equation {n — 1) we deduce that the only possible value for 
(3 is P = 0. Hence, we have the following result. 

Theorem 6 Let C be a real polynomial curve and let z(t) = c n t n + c n _i* n ~ 1 + 
• • • + cxt + Co be a polynomial, proper, real parametrization of C in complex 
form. If c n _i = 0, then the following statements hold: 

(i) If n = 3 (i.e. C is a cubic) then C is symmetric with respect to the point c . 
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(ii) If n > 4 then C has central symmetry if and only if c k = when k G 
{1, . . . , n — 2} is even. Moreover, in that case the symmetry center is Cq. 



Proof. In the case of (i), one can check that a — — 1, /? = fulfill the equations 
(1), (2), (3) of the system S. Then plugging a = —1,(3 = into (0) we 
get zq = Co. In the case of (ii), when we substitute (3 = in the equations 
(1), . . . , (n — 2) we get c k ■ (— l) fe + c k = for k = 1, . . . , n — 2. When k is 
odd this equality is clearly fulfilled, and when k is even the equality holds iff 
c fc = 0. □ 



To end this section, we provide information, in the following table, on some 
of the examples that we have tried. One may observe from the timings given 
that we can check curves of high degrees, algebraic coefficients, serious norms, 
etc. in seconds. The degrees and norms that we spell below correspond to ip(t) 
(not to the implicit equation, costly or even impossible to compute in some 
of these cases); furthermore, for each example we make explicit whether the 
curve shows central symmetry or not. The column of "norm" corresponds to 
the maximum of the infinity norms of the components of tp(t). In all the cases 
we show the timing (in seconds) corresponding to the method suggested by 
Theorem [5j to provide evidence of its efficiency; nevertheless, in some exam- 
ples, one can get a quicker response by detecting that (3 is not real (in which 
case this type of symmetry cannot exist) or by applying Theorem [6] The first 
column in the table below corresponds to the number of the example in our 
database. 



Ex. 


Deg. 


Norm 


Cent. Sym. 


Timing 


Comments 


2 


83 


2 81 


Yes 


0.718 












3 


3 


15 


Yes 


0.343 


Cubic with C2 = (see Th. 


6 


) 






4 


7 


15 


Yes 


0.437 


Curve of degree > 4 with c n _i = 


'see Th. 


6 


) 


37 


21 


2 30 


No 


0.468 


(3 £ R: with this, 0.172 sees. 








38 


45 


194 


No 


0.421 


(3 E: with this, 0.109 sees. 








39 


35 


8.06 • 10 9 


Yes 


1.248 


Algebraic coefficients. 










40 


77 


83 


No 


0.500 


(3 = but (—1,0) is not a solution of the syst 


em. 


41 


95 


5.89 • 10 27 


Yes 


5.850 


Algebraic coefficients. 
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3 Mirror Symmetry 



3. 1 Detecting Mirror Symmetry 

Along this section, as in the preceding one, we consider a real polynomial 
curve C furnished with a proper parametrization written in complex form as 
z{t) = x(t)+i-y(t). Also, we write the conjugate of z{t) as z(t) = x(t)—i-y(t). 
Now we say that C has mirror symmetry if there exists an axis C (called the 
symmetry axis) such that C is symmetric with respect to C. One can see 
that C has this type of symmetry if and only if there exists a movement 
M, composition of a translation and a rotation, such that M(C) is symmet- 
ric with respect to the x-axis. Hence, C has mirror symmetry if and only 
there exist Zq G C, and G [0,27r), such that z(t) = (z(t) — Zq) • e 1 '* is, in 
complex form, the parametrization of a curve C, symmetric with respect to 
the x-axis. Here, Zq defines the translation, and <j) defines the angle of rota- 
tion (see Figure 2); so, z(t) is a parametrization of M{C). Notice that if £ 
is a symmetry axis of C, then any point of £ can be chosen to be Zq . Now 
since the symmetric of a complex number a G C with respect to the x-axis 
is the conjugate a, it follows that C exhibits mirror symmetry if and only if 
w(t) = (z(t) — Zq) ■ e 1 '^ and w(t) = (z(t) — z ) ■ e 1 ^ both parametrize the same 
curve C. Furthermore, if z(t) corresponds to a proper parametrization then 
w(t) and w(t) define also proper parametrizations of the corresponding curves 
(because these parametrizations are obtained by composing a bijective pla- 
nar transformation with a proper parametrization). Taking also into account 
Lemma |3j the following theorem, analogous to Theorem [4] in Section [2j holds. 

t y 




X 



Fig. 2. Mirror Symmetry 

Theorem 7 The curve C exhibits mirror symmetry if and only if there exist 
a line C and a, (3 G R, a ^ 0, such that for any z G C, it holds that w(t) = 
w(at + (3). 
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Let us see how to take advantage of this theorem for efficiently checking, and 
compute in the affirmative case, the existence of mirror symmetry. We need 
first the following lemma. 

Lemma 8 A polynomial curve cannot have more than one symmetry axis. 

Proof. If C has m symmetry axes then by symmetry, each infinite branch 
gives us m infinite branches more. However, since C is polynomial then it has 
only two infinite branches, and therefore m — 1. n 

We will also need the following result, where we use the notation and termi- 
nology of Theorem [7j 

Lemma 9 Ifw(t) = w(at + (3) for a, (3 £R, then a = — 1. 

Proof. If w(t) = w(at + (3) for a, (3 E R, then by Theorem [7] the curve C 
is symmetric with respect to a certain axis, and w(t) parametrizes a curve C 
which is symmetric with respect to the x-axis. Let us write w(t) = (a p t p + 
a p _it p ~ 1 + - ■ ■ , b q t q + b q _ 1 t q ~ 1 + ■ ■ ■). Because of the symmetry of C with respect 
to the x-axis, we have that q must be odd (so that the sign of y(t) when t goes 
to oo or — oo is different). Now since w(at + (3) = w(t), then 

b g a q t q + ■■■ = -h q t q + ■■■ 



Hence, a q = —1. Since q is odd, b q ^ and a is real, then we deduce that 
a = —1. □ 

Now taking into account Theorem [7] and the expression of w(t), after some 
calculations we obtain the following system W, similar to the system S in 
Section H 

(c - zo) ■ e-^ = (o - zo) + ct/3 + c 2 /3 2 + • • ■ + e^/?"- 1 + c n (3 n 

c x ■ e -1 ' 2 ^ = da, + c 2 • 2a/3 H h c n _i(n - l)af3 n ~ 2 + c n na/3 n_1 

C 2 -C-^= C 2 a 2 +■■■+ C n _! ("- 1 > 2 ("- 2 > a 2 /3 n-3 + Cn r^ a 2pn-2 

r ■ p-' 1 ' 2 ^ — r a n 



We denote the above equations as (0), (1), ... , (n), respectively (i.e. the equa- 
tion (fc) is the one containing the power a k ). Notice that the equation (0) is 
the only one containing z$, and is linear in zq. Furthermore, from Theorem [7] 
one can see that the existence of mirror symmetry is equivalent to the exis- 
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tence of G [0, 2ir), a, f3 G K such that: (i) these values fulfill the equations 
(1), . . . , (n); (ii) for these values, the equation (0) corresponds to a real line 
(which would be the symmetry axis C). Now by Lemma [9j if W is consistent 
then a = — 1; also, by combining the last two equations we get that 



d = £(a) 



nlcj 2 



And since a = — 1, we get 

CnCn— 1 "I - C n C n _i 



0=e(-i; 



n C„ 



In the numerator of this expression we recognize c n c n _i + c n c n _i (i.e. the sum 
of a complex number and its conjugate), which is a real number; so, G R. 
Now in order to verify whether the equation (0) corresponds to a real line, the 
following results are useful. The first one can be proven in a straightforward 
way. 

Lemma 10 Any real line Ax + By + C = is transformed by means of the 

complex change x = . y = into + + C = 0, where 7 = 4 + i^-. 

2 _ _2i z a 

Conversely, any equality 72; + 72; + C = with C real corresponds to a real 

line. 



The above lemma is used for proving the following result. Here, we denote 
Qifi) = Cl /3 + c 2 f3 2 + ■■■ + c n /3 n , and 

1 1 

Proposition 11 The equation (0) of the system W corresponds to a real line 
if and only if Q*(/3) 6 1. 

Proof. After substituting — e~ l2<t> = e l ^~ 2 ^ = e l2 ^~^ in the equation (0), 
and dividing the whole equation by e 1 ^ - ^, we get 

e-W-Vzo + e^-^z = c ■ e^'® + c • + • e"^ 



is e 



i(7T/2-0) 



Now from Lemma 10, and taking into account that the conjugate of e 



-i(7T/2- 



it holds that the above expression corresponds to a real line iff the 
right hand-side of the equation is a real number. However, one may see that 



c • e 



»(f- 



+ c • e 



-Ki- 



ll 



is already real, because it is the sum of a complex number, and its conjugate. 
So, we get a real line iff Q(j3) ■ e" 1 ^"^) is a real number. Finally, using that 
_ e -i2<j> _ e i2(|-</>) an( j e q Ua tion (n), one has that 



e -i(f-0 = 1 = 1 = 1 = ./ ~ c " 

e i(f-0) y/- e -i2<t> V^l- 1 )" 



Multiplying and dividing by c n into the square root, the result follows easily. 

□ 

So, we can put together all these observations to characterize the existence of 
mirror symmetry. For this purpose, it is preferable to consider a new system 
W, which is obtained from W by substituting 

Cn 

(from the last equation) into the n first equations, and dividing out each 
equation (k) by a k : 

C n ■ [(C - Z ) + Ci/3 + C 2 /3 2 + • • • + Cn.^- 1 + C n fi n ] 

c n -[d + 2c 2 (3 +•••+ c n . 1 (n-l)f3 n ~ 2 + 

Cn ■ [C2 + ■ ■ ■ + Cn-l ' 1 " 1 ;""^ " 3 + Cn^^ 2 ] 
C n _iC n • a = Cn - [c n -l + C n nf5] 



(c - z )c n ■ a n = 
c\c n ■ a n ~ l = 
c 2 c n ■ a n ~ 2 = 



We denote these equations as [0], [1], . . . , [n — 1], respectively. Furthermore, 
we write the equation [0] as r(z ,z ,f3) = 0. Then the following result, which 
can be deduced from the results in this section, holds. Observe that this result 
provides a characterization for the existence of mirror symmetry, and a closed 
expression for the symmetry axis, when it exists. 

Theorem 12 The curve C has mirror symmetry if and only if the following 
two conditions hold: (1) a = —1, (3 = £(— 1) fulfill the equations [1], . . . , [n — 1]; 
(2) Q*(f3) € K. Moreover, in that case the symmetry axis of C (in complex 
form) is r(z,z, (3), i.e. 

c n z - c n (-l) n z + c c n (-l) n - c n c - (ci/3 + c 2 /3 2 + ■■■ + c n (3 n ) = 

Remark 2 One can check that the conditions (1) and (2) are independent, in 
general. Take for instance a cubic parametrized by z(t) = t 3 + c 2 t 2 + c±t + c , 
i.e. c 3 = 1. In this case, one can check that a = — 1 and (3 = £(— 1) = 
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— o( c 2 + c 2 ) IM-fitt equations [1], [2] iff Ci = c\ — -c\ + -c\. However, 

o o o 

Q*(/3) = - — [ci + c 2 ) • (9ci - 2cj - c 2 c 2 + c\). So, for example c x — 1 + 2/3i, 

c 2 = C3 = 1 satisfy [1], [2] but the value of Q*{(3) in that case is —lAi/27, 
which is not real. 

Remark 3 One may see that if C has a symmetry axis C, then both must 
intersect. Indeed, since C is symmetric with respect to C, if it has a point in 
the half-plane over C then it must have also another point in the half-plane 
below C But since C is a continuous curve (i.e. there is no t-value where it 
becomes infinite), both points must be connected, and therefore the symmetry 
axis must be crossed. 



Theorem 12 provides also the following result on the existence of symmetries 
ofC. 

Proposition 13 The following statements are true: 

(1) If c n — (— l) n c n = ; and C has a symmetry axis, then it is parallel to the 
x-axis. 

(2) If (— l) n c n + c n = 0, and C has a symmetry axis, then it is parallel to the 
y-axis. 

Proof. Let us see (1). For this purpose, assume that £ is a symmetry axis for 
C. The intersection of L with the x-axis can be found by substituting z = z 



in the closed expression for z provided in the statement of Theorem 12 This 
substitution yields (c n — (— l) n c n )z + • • • = (where • • • comprises terms 
where z is not present), and hence when c n — (— l) n c n = no value for z can 
be found. So, (1) follows. The statement (2) is proven in a similar way but 
imposing z = —z , instead of z = z . □ 

The above proposition can be used for proving the following, clearer, result on 
the existence of symmetries. Here we denote deg t (x(t)) = r, deg t (y(t)) = s. 

Theorem 14 The following statements hold: 

(i) If r > s and r is even, and C has a symmetry axis, then it is parallel to the 
x-axis. 

(ii) If r > s and r is odd, and C has a symmetry axis, then it is parallel to the 
y-axis. 

(Hi) If r < s and s is odd, and C has a symmetry axis, then it is parallel to the 
x-axis. 

(iv) If r < s and s is even, and C has a symmetry axis, then it is parallel to the 
y-axis. 
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Proof. If r > s then n 



and c„ is real, i.e. c„ = c„. Now if n, i.e. r, is 



even we are in the hypotheses of statement (i), and the result follows from the 



statement (1) of Proposition 13 ; if n, i.e. r, is odd we we are in the hypotheses of 



statement (ii), which follows from the statement (2) of Proposition 13 Finally 
if r < s then n — s, and c n is a purely imaginary number, i.e. c n = —c n ; then 
we argue as before. □ 



Theorem 12 leads to an algorithm for computing the symmetry axis of C (if 
any). This algorithm is illustrated in the following example. 

Example 2. Consider the curve C parametrized by (x(t),y(t)) , where 



x(t) = (2t + l) 20 + (2t + l) 18 + (2t + l) 10 + 1 + (2t + l) 21 - 3(2t + l) 5 + (2t + l) 3 
y(t) = -{2t + l) 20 - (2t + l) 18 - (2t + l) 10 - 1 + (2t + l) 21 - 3(2t + l) 5 + (2t + l) 3 



27ms curve is proper, and has been constructed starting from a simpler curve, 
then applying a change of parameters and finally rotating the curve j radians. 
One can check that the implicit equation f(x,y) of C has degree 21, that the 
infinity norm of this implicit equation is close to 2 500 , and that the absolute 
values of the leading coefficients of f w.r.t. x andy are equal and of size 2 441 . 
Using our method, it takes 1.388 seconds to construct the system W and check 
that: 



a 



■ 1, (3 = fulfill the equations [1],.. . , [n-1], andQ*((3) = 2^2 e 



R. So, from Theorem 12, C has one symmetry axis. 
• The equation of the symmetry axis is obtained by substituting the above 
values for a, (3 in the equation [0], which yields iz — z = 0, i.e. y = x. 

As in the preceding section, in the = we can provide a sharper 

result. In this case (3 = and therefore Q*(f3) = 0; so, the second condition in 
Theorem [12] always holds. In addition to this, the equations of W have the 
form mfc(a) = c n Ck<y k — c n Ck = for k = 1, . . . , n — 2. So, when a = — 1 we 
have 



m k (a 



c n Ck — c n Ck if k is even 
-(c„c fe + c n c k ) if k is odd 



Hence, the next result, which essentially follows from the above expression, 
follows. 

Theorem 15 Assume that c ra _i = 0. Then C presents mirror symmetry if and 
only if one of the following two conditions hold: (i) for every k = 1, . . . , n — 2, 
c n ■ Ck G M; (ii) for every k G {1, . . . , n — 2}, c n ■ c k G K when k is even, and 
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c n • Ck is a pure imaginary number (i.e. with null real part) when k is odd. 
3.2 Some More Prohibitions and Observations 

If C exhibits mirror symmetry and £ is a symmetry axis, then almost all 
elements in the family of lines Ax + By + C = normal to C must intersect 
C at an even number of points (the finitely many exceptions will be those 
lines that contain some point of C D C). We can exploit this for deducing the 
non-existence of mirror symmetry in some cases. In order to do this, we need 
the following previous result. 

Lemma 16 Let p(t) , q(t) be two real polynomials, and let A,B,C G R. Also, 
let RA,B,c(t) — Ap(t) + Bq(t) + C, and let R'abc(^) denote its derivative with 
respect to t. Then there does not exist (A ,B ) G K 2 such that RA ,B ,c(t) has 
multiple roots for almost all C 6 1. 

Proof. For a particular value of C G R, Ra ,b ,c{^) has multiple roots iff 
RA ,B ,c(t) = A p(t) + B q(t) + C and R' AoyB{hC (t) = A p'(t) + B q'{t) have a 
common factor. Moreover, R Ao , Bo ,c{t) has multiple roots for almost all C G K 
iff the former polynomials have a common factor regardless of the value of C. 
But this cannot happen because the second polynomial does not depend on 
C. □ 



Then we have the following result. 

Theorem 17 Let r = deg t (x(t)), s 
hold: 



= deg t (y(t)). The following statements 



(1) If r (resp. s) is odd, then C cannot have any symmetry axis parallel to the 
x-axis (resp. y-axis). 

(2) If r = s and r, s are odd, then C has at most one symmetry axis Ax + By + 
C = parallel to the vector (a r ,a s ), where a r ,a s are the leading coefficients 
of x(t) ,y{t), respectively. 

Proof. Let us see (1). If C has a symmetry axis parallel to the x-axis (resp. 
the y-axis), then the family of lines x = a (resp. y — b) must intersect C at an 
even number of points, for almost all a G ffi (resp. b G R); however, if r (resp. 
s) is odd this cannot happen. Now in order to see (2), one observes that if 
C has a symmetry axis C, then the family of normal lines to C has the form 
A x + B y + C = for a certain (A , B ) G M 2 , where C is a real parameter. 
The intersection of this family with C amounts to Ax(t) + By(t) + C = 0, and 
this equation must have an even number of solutions for almost all CGK. Now 
assume that the degree of this equation is odd, and let S stand for the number 
of real solutions of the equation; in that case, either A a r + B a s = 0, or 
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the equation has multiple solutions. However, by Lemma 16 the latter cannot 
happen for almost all values of C, and (2) follows. n 



Example 2, in the subsection 2, provides an example for the statement (2) of 
this theorem. Now by putting together Theorem 14 and Theorem 17 we get 



the following corollary of Theorem 17 



Corollary 18 If both r = deg t (x(t)) and s = deg t (y(t)) are odd, and r ^ s, 
then C does not exhibit mirror symmetry. 



The results in Theorem 14, Theorem 17 and Corollary 18 are displayed in the 



table below, where the notation of Theorem [17] is maintained. Observe that 
the results in the preceding sections do not say anything for the case when 
both x(t) and y(t) have equal and even degrees. 





r < s : 


No mirror symmetry 


r odd, s odd 


r = s : 


Axis parallel to (a r ,a s ), if any. 




r > s 


No mirror symmetry 


r odd, s even 


Axis parallel to y-axis, if any. 


r even, s odd 


Axis parallel to x-axis, if any. 




r < s : 


Axis parallel to y-axis, if any. 


r even, s even 


r = s : 


(nothing to say) 




r > s : 


Axis parallel to x-axis, if any. 



Finally, one might wonder whether central symmetry and mirror symmetry 
can coincide at the same time. Let us see that the answer, in our case, is 
negative. For this purpose, we first need the following result. 



Proposition 19 Let C be an algebraic curve, and assume that it presents 
rotation symmetry, with center P and angle 8. If £ is an axis of symmetry of 
C, then Pq G C 

Proof. Assume by contradiction that Pq C, and let Pq ^ Pq be the sym- 
metric of Pq w.r.t. C. Furthermore, let C be the symmetric of C w.r.t. C. Since 
C is by hypothesis symmetric w.r.t. C, we have C = C. Hence, a rotation of C 
around Pq with angle 9 amounts to a rotation of C around Pq with angle —9. 
But since C = C, it holds that C has two symmetry centers, namely Pq and 
Pq. However this cannot happen because C is algebraic, and therefore it can 
have just one symmetry center (see Theorem 5.3 in [B]). So, Pq G C. □ 



Now we can state the result. 
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Theorem 20 Polynomial curves cannot exhibit central symmetry and mirror 
symmetry at the same time. 



Proof. If C has a center of symmetry P an d a symmetry axis £, then by 
Proposition 19 it holds that P e C. In that case, it is easy to see that C must 
also have another symmetry axis £', namely the normal line to £ at Pq. But 
by Lemma [HI this is not possible. □ 



3.3 More Examples 



We finish this section with a table showing data on some of the examples that 
we have tried. From this table, and as it also happens in the case of central 
symmetry, we can see that the method provided here can deal efficiently with 
serious inputs. In all the cases we provide the timing (in sees.) of the algorithm 



deriving from Theorem 12 however, under certain circumstances (i.e. when 
Q*{0) ^ K, when c n _i = 0, or in some cases corresponding to the table in 



Subsection 3.2) we can use some of the results in the preceding subsections to 



derive an answer faster. 



Ex. 



Degree 



Norm 



Mirr. Sym. 



Timing 



Comments 



6 
7 
8 
9 

10 
11 
12 

42 
43 
44 



21 
21 
69 
45 
45 
60 
91 
77 
35 
56 



2 31 

3.5 
276 
461 
194 
844 

2 89 

2 71 

2 38 
2 52 



Yes 
Yes 
No 
No 
No 
Yes 
Yes 
No 
No 
Yes 



1.388 
0.983 
3.229 
1.653 
1.653 
2.465 
15.803 
2.433 
2.901 
23.385 



Case c n —i 
Case c n _i = (see Th. 



Algebraic coefficients. 
(see Th. 



15 



15 



with this, 0.203 sees.) 
with this, 0.203 sees.) 



r even, s odd (see Th. 



14) 



Q*(f3) £ K; with this, 2.293 sees, 
r, s odd (see Th. 
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with this, sees. 



Algebraic coefficients. 



4 The Floating Point Case 



In the preceding sections, we were assuming that the input was given exactly. 
However, it can happen that the coefficients are known only up to a certain 
precision e, that can be understood as the number of correct decimal positions 
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of the coefficients. We discuss this situation here. The overall idea is that while 
in the exact case one first tests whether there exists a certain type of symmetry, 
and one computes the symmetry center or the symmetry axis when the answer 
is affirmative, here it is better to proceed just the converse way: one uses the 
results in the sections before to provide the temptative symmetry center, or 
symmetry axis, and afterwards one evaluates whether the curve behaves in 
approximately a symmetric way with respect to the computed point/line. 



4-1 Central Symmetry 



One can use the closed formula provided by Theorem [5] also in this case. How- 
ever, the inaccuracy of the coefficients causes the estimation of the symmetry 
center by this formula to be also inaccurate. In order to measure this inac- 
curacy, we can differentiate the formula for zq provided in Theorem [5j and 
interpret the norm of the differential as a measure for the error the value z 
is known with. Hence, we have that 

dz = dc + \{dc x p + dc 2 /3 2 + ■■■ + dc n /3 n ) + \{c x + 2c 2 /3 + ■■■ + nc n (3 n - r )dl3 

_ _ 



dc n — x - On C n _ \ 'dCn 



Now, let |c| 
for all i = 0, 1, . 



Furthermore, d/3 = '■ 
of c; then we have that |q| < 
i = 1, . . . , n. Hence, substituting the expression for d(3 in the above expression, 
taking norms, and making use of usual properties on norms, we have that 



denote the infinity norm 
. . , n. Also, \dci\ < e for 



\dzo\ < e 



I + [-n + l)max{l,\/3\ n } 



Notice here that \(3\ = ( ^ Cw ~^ Cn ^ , and hence for a fixed value of \c n -i/c n \ 
the above bound grows linearly as n goes to oo. This bound, however, has two 
drawbacks: (i) it can be small although central symmetry is not present; (ii) 
it can be slightly high even when central symmetry is present. So, in practice 
it is only useful when it is small, and one knows that the curve has this kind 
of symmetry. An alternative possibility for evaluating central symmetry is to 
take into account the test provided by Theorem |5j i.e. checking whether the 
equations (1), ...,(&) of the system S (see Subsection 2.2) are fulfilled. These 
equations correspond to polynomial functions fk(c ,Ci, . . . ,c n ,a, (3), that in 



presence of symmetry should vanish when evaluated at a = — 1,(3 — 
and the coefficients of the parametrization. In our case, because of the limited 
precision the q's are known with, even in presence of approximate central 
symmetry these functions do not evaluate as 0. However, experimentation 
shows (see later) that when the curve exhibits central symmetry these values 
tend to be small (whenever the parametrization is not too big, in norm or 
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degree). One can prove (by proceeding as before) an upper bound 0{e ■ 2 n ■ n ■ 
max{l, |/3| n }) for the inaccuracy of these evaluations, but in our experiments 
we have observed that this bound is generally too high. A measure like L = 
max{|/i|, . . . , \ f n \}, where the JVs stand for the evaluations of the /fc's, gives 
a better idea. Another possibility is to evaluate the central symmetry with 
respect to zq in a statistical fashion. For this purpose, we observe that from 



the results of Subsection 2.2 , the symmetric of <p(t) with respect to the point Zq 
(if this is really a symmetry center) is <p(—t + /3). So, we generate 400 t-values, 
giving rise to points P t = <^(i)'s, and we compute the P t = <p(—t + /3)'s; then 
for each t- value we compute the difference between the distance d(P t , Pt), and 
the value of 2 • d(P t ,Za) (i.e. twice the distance d(P t ,zo)). We represent this 
difference by 8P t , and its relative value (expressed as a percentage) by SP t , 
i.e. 5P t = ^(pfzo) ' -'■OO. Finally, the arithmetic mean M (after taking out the 
outliers) of these values measures whether z is really a symmetry center, or 
not. 



In the table below we provide some information concerning 10 of the examples 
that we have tried. In each case, we spell: the number of the example, the 
degree of the parametrization (i.e. the maximum power of t appearing in it), 
the infinity norm of the parametrization, the precision (i.e. the accuracy of 
the coefficients), and the values of the three above parameters: |g?Zo|j L, M. 
Furthermore, in each case we have explored the curve graphically and we have 
confirmed the existence or not of approximate symmetry, 



Ex. 


Degree 


Norm 


Prec. 


Symm. 


\dz \ 


L 


M 


16 


13 


1964.14 


10" 6 


Yes 


0.033 


0.001 


0.0026 


17 


5 


28.41 


10~ 3 


Yes 


0.273 


2.23 • 10~ 7 


3.35 • 10~ 7 


18 


9 


3.0006 


10" 4 


No 


0.015 


42.82 


9.18 


19 


29 


3.0006 


10~ 4 


No 


0.0045 


14.85 


9.70 


20 


29 


10.002 


10" 4 


No 


0.0045 


26.742 


12.94 


21 


25 


3.002 


10" 4 


Yes 


0.0039 


0. 


5.40 ■ IO" 8 


22 


39 


14.0028 


10" 4 


No 


0.0061 


31.25 


21.47 


23 


57 


182.73 


10" 4 


Yes 


0.00875 


0.0076 


1.75 • 10" 4 


24a 


53 


2 50 


10~ 6 


Yes 


6.11 • 10 12 


1.8 • 10 11 


4.11 • 10 15 


24b 


53 


2 50 


io- 16 


Yes 


611.28 


7.44 • 10~ 15 


10" 4 



It is interesting to observe that in examples like 16 or 24b, the bound for \dzo\ 
works worse than L or M. On the other hand, one may notice several examples 
(18, 19, 20, 22) where \dz \ is small, although there is not central symmetry. 
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Also, the example 24a shows a curve with very high coefficients, that we have 
approximated with precision 10~ 6 ; even though the curve is clearly symmetric 
(it is a small perturbation of a symmetric curve with algebraic coefficients), 
the three parameters were very big and the symmetry was not recognized by 
any of them. We needed to use much higher precision (10 -16 ) to recognize the 
symmetry (see Example 24b). 



4-2 Mirror Symmetry 



In the exact case, one computes the symmetry axis C by using the (0) equation 
of the system W'. However, in the approximate case C is better determined 
by separately computing 0, and a point P of £, written as z in complex 
form. Unless the axis is parallel to the x-axis (or c n — c n (— l) n is very small 
in norm) we choose Pq to be the intersection of C with the x-axis (in complex 
form, this will correspond to a real number, therefore satisfying z = z ); if 
£ is parallel to the x-axis (or c n — c n (— l) n is small), we will compute instead 
the intersection with the y-axis (a purely imaginary number, hence fulfilling 
zq = —zq). In order to compute 0, and whenever we are in a case where 
the value of cannot be determined from the beginning (see the table in 



Subsection 3.2), we will use the equation [n] of W; so, we can easily see that 



where r n = c ^ ' . This computation is well-conditioned whenever Re(r n ) is 
not too small. Afterwards, we can compute Zq from the equation [0]: here we 
will get a formula for z with c n + 5(—l) n c n in the denominator, with 5 = — 1 
whenever we impose z = z , and 5 = 1 whenever zq = — z . As a consequence, 
if this denominator is small in norm we have a bad-conditioned case. A general 
analysis of the sensitivity of 0, zq in this case is more complicated than in the 
case of central symmetry, and involves not only general parameters (precision, 
degree, norm, etc.) but also other questions, like the closeness of \c n +5(— l) n c n \ 
to 0. So, in this case we suggest two possible ways for assessing the existence 
of symmetry with respect to £: 

(i) The equations of the polynomial W correspond to polynomial expressions 
9k(co, Ci, . . . , c n , a, (3) = 0. Hence, a first indicator of the existence or not of 
symmetry is the maximum absolute value L of the g^s, when evaluated at 
a = —1, (3 = £(— 1), and the parametrization. 

(ii) Statistical criteria. In our experiments we have used the following: we gen- 
erate points of C (by substituting 400 values of t in z{t)). If C is really a 
symmetry axis, one can check that the symmetric of z(t) with respect to C 
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is 



z(t) = e" 2 ^ • (z(t) - zo) + z = . (z(t) - zo) + z 

Thus, denoting the distance of z(t) to C as d t , and the distance between 

z(t) and z(t) as dt, one might consider the arithmetic mean of the relative 

\dt — 2cm . . . 

errors o t = ~ , m percentage. In our case, we improved this a little 

d t 

by: (a) computing the intersection points t\, . . . , t r of C and C; (b) giving 
400 values around each of these points (i.e. 400 for the first one, another 
400 for the second, etc.); (c) computing the arithmetic means X\,...,x r of 
the deviations of the points obtained around each tf, (d) finally, taking the 
median of x±, . . . ,x r . We decided to compute the intersections in order to 
avoid using points with very high distances from C, which might lead to 
numerical problems. Also, we took the median of X\, . . . , x r to avoid errors 
due to fake intersections detected by the Maple solver, which we detected 
specially for high degrees/norms. The statistical parameter produced this 
way is denoted, in the following, as M. 

The problem in the case of (i) is that L can be big even when the curve is 
approximately symmetric, due to inaccuracies. In the case of M, we found 
that some polynomial curves which are non symmetric when examined "lo- 
caly" , show however small (and even very small) values for M. This suggests 
a "global" symmetry, i.e. symmetry from a statistical point of view, which is 
what M is really evaluating. We see an example of this in Figure 3. At the 
left, we have the curve 

(p(t) = (-73.073 + 97.097t - 62.062t 3 - 56.056t 9 + 87.087t 10 + l.OOlt 14 + LOOK 15 , 
-17.017 + 71.07H 5 - 44.044t 6 + 80.080t 8 - 82.082t n + l.OOlt 14 + 62.062t 15 ) 

plotted for t G [—1,1], and at the right we have the same curve for t G 
[—10,10]. We see that the curve is not "locally" symmetric, although "glob- 
ally" it is. For this curve, we have L = 194 and M = 5.973 • 10~ 5 , therefore 
showing global symmetry. 

In the following table we provide information on some of the examples that 
we have tried. We see that only two of the 10 curves shown in the table show a 
high value of M; so, we must conclude a common tendency between polynomial 
curves to mirror symmetry. Also, we find several curves with big values for 
L, but small values for M. Finally, it is worth observing that the examples 
25, 29, 30, 34, 37 correspond to symmetric curves that have been slightly 
perturbed, while the other examples in the table correspond to perturbations 
of non-symmetric curves. Hence, we can see that in general the closeness to 
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Fig. 3. Global Symmetry 

"local" symmetry is well recognized; an exception is the curve 34, because of 
its high coefficients (check the norm column for this curve). 



Example 


Degree 


Norm 


Precision 


L 


M 


25 


21 


4.5 


10 -io 


lO" 9 


0.00003 


26 


69 


275.99 


10" 4 


663.16 


52.53 


27 


45 


461 


10" 4 


970 


0.01087 


28 


8 


97 


10" 4 


192.51 


270.576 


29 


60 


843.991 


0.001 


0.002 


0.0 


30 


30 


3.003 


0.01 


3.27 


0.01159 


31 


51 


5.005 


0.01 


11.55 


0.25 


34 


21 


2.59 ■ 10 9 


0.01 


17588.46531 


0.00018 


35 


6 


120.32 


0.01 


266.45 


4.03558 


36 


17 


97.097 


0.01 


166.166 


2.64- HT 8 


37 


20 


3.006 


0.01 


3.673 


2.95707 



5 Conclusions and Further Work 

Here we have presented an efficient method for detecting the symmetries, and 
computing them in the affirmative case, of algebraic curves defined by means of 
polynomial parametrizations. The method is a combination of geometric and 
algebraic ideas: on the one hand, we observe that the nature of the symmetry 
leads to a new parametrization of the curve; on the other hand, whenever we 
start from a proper parametrization, this second parametrization must also 
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be proper, and therefore is related with the first one by means of a certain 
algebraic relationship (see Lemma [3J. This way we obtain a polynomial sys- 
tem with three unknowns but with a triangular structure, that can therefore 
be solved in a fast and efficient way. In fact, we provide closed expressions 
for the temptative solutions of the system, and for the elements defining the 
symmetries of C, if any. Furthermore, we have also reported the experimen- 
tation done in the case of floating point inputs, showing that our results can 
be also useful in that case. Now it is natural to wonder if this method can be 
extended to other situations. In the case of rational, not simply polynomial, 
parametrizations of plane curves, we can apply the same strategy to build 
proper parametrizations, which however are related by means of a linear ra- 
tional change of coordinates (see Lemma 4.17 in [9]), and not simply a linear 
one; this leads to a polynomial system with more unknowns than in the poly- 
nomial case, where the triangular structure is lost. In the case of space curves, 
the difficulty comes from the fact that space rotations are harder to treat; in 
particular complex numbers cannot be used anymore, although they can be 
replaced by quaternions (see for example j3]). Finally, in the case of surfaces 
there are also results relating proper parametrizations of a same surface (see 
[8]); however, here the number of parameters increases, which leads to a more 
difficult framework. So, in all these cases the method must be adapted. Our 
current efforts are focused to adapt the ideas presented in this paper to these 
situations, which can therefore be considered as on-going work. 
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